Unitarity plus causality implies localizability 
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We consider a graph with a single quantum system at each node. The entire compound system 
evolves in discrete time steps by iterating a global evolution U. We require that this global evolution 
U be unitary, in accordance with quantum theory, and that this global evolution U be causal, in 
accordance with special relativity. By causal we mean that information can only ever be transmitted 
at a bounded speed, the speed bound being quite naturally that of one edge of the underlying graph 
per iteration of U. We show that under these conditions the operator U can be implemented locally; 
i.e. it can be put into the form of a quantum circuit made up with more elementary operators — 
each acting solely upon neighbouring nodes. We take quantum cellular automata as an example 
application of this representation theorem: this analysis bridges the gap between the axiomatic and 
the constructive approaches to defining QCA. 



I. MOTIVATIONS 

A physical system is described in quantum theory by a 
state vectors \ip) (a unit vector in a Hilbert space TL), and 
evolves from time t to time t' according to a unitary op- 
erator U. The definition of a unitary operator is made in 
general terms: induce an bijection between two orthonor- 
mal bases (WU = WU = I). But what if we want a more 
hands-on, operational description of a unitary operator? 
In finite dimension we know that they can be spectrally 
decomposed (U = ^ x e, lX *\cf) x ){4> x \), but also that they 
can be approximated up to arbitrary precision by a cir- 
cuit composed of the universal quantum gates H, Phase, 
C — Not — via the Solovay-Kitaev theorem. However in 
infinite dimensions spectral theory becomes quite compli- 
cated, and nothing tells us whether the operator can be 
expressed as a quantum circuit. Often it can be difficult 
to provide an operational description of unitary opera- 
tors over infinite dimensional spaces; in that sense they 
remain abstract mathematical objects. 
Usually the infinite number of degrees of freedom arises 
from a position degree of freedom, i.e. when space comes 
into the picture. The canonical example is that of the 
wave function of a particle on a line. But physics then 
tells us something else about the evolution, namely that 
if the particle is well-localized within a region R at time 
t', it was not to be found outside the region R ± c(t' — t) 
at time t. This is a case of causality, and what causality 
says in general is that if we distinguish different "places" , 
some of them close to one another, some of them distant, 
and if the interval (t' — t) is sufficiently small, then the 
state associated to some place x at time t' should only 
depend upon the state associated to the neighbours of x 
at time t. 
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A composite physical system AB is described in quan- 
tum theory by a state vector \tp) AB (a unit vector in 
H A <8> T~t B ), which in general cannot be decomposed into 
state vectors \ip) A /\ip) B associated to subsystems A/B — 
due of course to entanglement. In order to still be able 
to speak of "state associated to some place" we must 
switch to the well established formalism of states and 
partial traces. So now a composite physical system is de- 
scribed by a state p AB (i.e. a unit trace positive operator 
over Ti. A C8) Tt B ), and the state associated to place A is 
p\a = Trs(p). And hence we can express causality by 
saying that if p taken to p' over a short enough period 
of time, then p'\ x should be a function of p|a/^i with Af x 
designating the neighbours of x. 

In this paper we will study unitary causal operators for 
they own sake, and we will show that they have a lot of 
structure. For instance the inverse of a unitary causal 
operator is also a unitary causal operator. More impor- 
tantly, we will show that they decompose into a prod- 
uct of local operators, i.e. operators which act solely 
upon neighbourhoods. This structure theorem is a gen- 
eral representation theorem for unitary causal operators, 
which yields an operational description of them — just 
like the two aforementioned representation theorems did 
for finite-dimensional unitary operators. 

This structure theorem follows a tradition of postu- 
lating causality over a global dynamics and then show- 
ing that the dynamics can be implemented locally - 
which is the difficult direction to go to of course (mainly 
due to entanglement), the converse direction being al- 
ways trivial. It provides a general answer to this ques- 
tion under the unitarity and discreteness hypotheses, 
encompassing (up to some details) the previous results 

on this issue by Beckman et Al. j3[ for two systems, 

Schumacher and Westmoreland [2J| for three systems, 
and Schumacher and Werner for a line of translation- 
invariant one-dimensional systems (23j , i.e. Quantum 
cellular automata (QCA). As regards QCA more specif- 
ically, a corollary of our results is that the axiomatic 
definition of ?vdimensional QCA proposed in p3j does 
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admit an operational. Block-structured QCA description 
after all — and conversely that these seemingly specific 
Block-structured QCA |9|] are in fact general instances 
axiomatic definition. 

We start with basic definitions and properties (Sections 
HTl and before we prove our main result in Section 
IIV1 We then discuss related works mainly in quantum 
cellular automata, where we take our inspirations from, 
and discuss the more general perspectives (Sections IVl 
and [Vil- 
li. DEFINITIONS 

So far in this intuitive motivation towards the concept 
of unitary causal operator we have been speaking about 
"places, some of them close to one another, some of them 
distant" . As we seek to capture this idea in the most gen- 
eral, and yet simple and formal manner, we shall identify 
those "places" with the nodes of an arbitrary graph, and 
say that two nodes are "close" whenever they are related 
by an edge. We first need to make rigorous the idea of a 
graph, with a quantum system at each node. 

Definition 1 (quantum labeled graph) 

A quantum labeled graph (QLG) is a tuple T = (V,£,TL) 
with: 

- V, the nodes (a countable set); 

- £, the edges (a subset ofVx V); 

- TL, the labels (a countable set of Hilbert spaces). 

We denote by Af x = {y \ (x, y) G £} the set of direct 
neighbours of the node x, with x an integer ranging over 
V. 

To each node x there is an associated alphabet Yj x and 
hence an Hilbert space TL X = TLs* . Now say the graph is 
infinite. The difficulty here is that as we have mentioned 
an infinite tensor product Hilbert spaces "® N TL X " is in 
general not a Hilbert space, so we must take the following 
detour: 

Definition 2 ((finite) configurations) 

A (finite) configuration c of a QLG T = (V,£,TL) is a 
function c : N — ► N, with x i — ► c(x) = c x , such that: 

- c x belongs to T, x ; 

- the set {x\c x ^ q} is finite. 

The set of all finite configurations of a QLG be denoted 
Cf again. 

The idea is that finite configurations are the basic states 
of the quantum systems labelling the graph. The follow- 
ing definition works because C/ is countable: 

Definition 3 (superpositions of configurations) 

We define TLc f be the Hilbert space of configurations of 
a QLG r = (y,£,H), as follows: to each finite config- 
urations c is associated a unit vector \c) , such that the 
family (|c)) ceC is an orthonormal basis ofTLc s - A state 
vector is a unit vector \?p} in TLc f - A state is a trace-one 
positive operator p over TLc f ■ 



Note that TLc f is entirely defined by the set of Hilbert 
spaces TL = (TL X ). From now on we will write TL instead of 
TLc f ■ Note also that the state p captures the state of the 
entire compound system, whereas p\ x stands for the state 
which labels node x of the graph, where we introduce the 
notation A\ s for the matrix Tr A n but the systems in s(A). 

Definition 4 (Causality) 

A linear operator U : TL — ► TL is said to be causal with 
respect to a quantum labeled graph T if and only if for 
any p, p' two states over TL, and for any x G Z, we have 

pW.=/W. ^u(p)u^\ x = u(p')u^\ x . (i) 

In other words: to know the state of node number x, 
we only need to know the neighbouring of nodes M x . 
Unitarity is as usual: 

Definition 5 (Unitarity) 

A linear operator U : TL — ► TL is unitary if and only if 
{U\c) | c G Cf} is an orthonormal basis of Tic r 

Hence we have defined the main object of our discourse: 
unitary causal operators. This concept of unitary causal 
operator generalizes the two-systems definition by Beck- 
man, Gottesman, Nielsen, Preskill Q and the three- 
systems definition by Schumacher and Westmoreland 



III. PROPERTIES 

Let us begin by proving some fundamental facts about 
unitary causal operators, which may also be regarded as 
alternative formulations of causality. Proposition [T] ex- 
presses causality in the Heisenberg picture, as a condition 
on the evolution of observables. Whenever we say that a 
linear operator A is localized upon a region R, we mean 
that A is of the form An ® Iy\i?, i.e. it is the identity 
over anything that lies outside of R. Morally, A is an 
observable in the following result. 

Proposition 1 (Dual causality) 

Let U be a causal linear operator with respect to a quan- 
tum labeled graph T. This is equivalent to saying that for 
every operator A localized upon node x, then U* AU is 
localized upon the nodes in Af x . 

Proof. [=>]. Suppose causality and let A be an oper- 
ator localized upon node x. For every states p and p' 
such that pl./^ = p'\u x > we have (UpW) \ x = (Up'W) \ x 
and hence Tr (AU pW) = Tr(AUp'W). We thus get 
Tr (WAUp) = Tr (TPAUp 1 ) . Since this equality holds 
for every p and p' such that p\tf m = p'|a4, what we are 
saying is that the W AU does not discriminate differ- 
ences between p and p' whenever they lie outside of N x . 
In other words U'AU is localized on the nodes in M x . 
[•£=]■ Suppose dual causality and p\j^ x — p'\j\T x - Then, 
for every operator B localized upon the nodes in Af x , 
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Tr (Bp) = Tr(Bp'), and so for every operator A local- 
ized upon node x, we get: Tr (AUpW) = Tr (WAUp) 
Tr (WAUp') = Tr (AUp'W) . This entails (UpW) \ x = 
(Up'W) \ x . □ Proposition [2] expresses causality in terms 
of the inverse of the unitary causal operator U. When- 
ever we speak about the transpose of a quantum labeled 
graph r, we mean as usual the quantum labeled graph 
r T which is obtained just by changing the direction of 
the edges. The neighbours of x in T T are designated by 
Mj. 

Proposition 2 (Inverse causality) 

Let U be a causal linear operator with respect to a quan- 
tum labeled graph T. Then W is a causal operator with 
respect to the transposed quantum labeled graph T T . 

Proof. Suppose causality, let A be an operator localized 
upon node x, and choose M an operator localized upon 
a node y which does not lie in Af x . That way x does not 
belong to Af y . But according to Proposition [1] we know 
that U'MU is localized upon Af y , and hence U'MU com- 
mutes with A. Now A i— > U AW is a morphism because 
AB i ► UAWUBW = UABW, and so via this mor- 
phism we can also say that UWMUW = M commutes 
with U AW . An since M can be chosen amongst to full 
matrix algebra Md(C) of the node y, this entails that 
UAW must be the identity upon this node. The same 
can be said of any node outside NJ '. So U AW is localized 
upon M x and we can conclude our proof via Proposition 

CD □ 

We will use both these propositions in order to establish 
our representation main theorem. 

IV. REPRESENTATION 

We will now show that unitary causal operators are 
implementable locally; i.e. that they can be put into the 
form of a quantum circuit made up with more elementary 
operators - each acting solely upon neighbouring nodes. 

Theorem 1 (Local representation) 

Let U be a unitary causal operator with respect to a 
quantum labeled graph T = (V,£,H). Then there exists 
D, (K x ), E, and \<f>) such that for all \ip), 

((g) D)([[ *r x )((g) em = |0) ® uty) 

where: 

• (K x ) is a collection of commuting unitary opera- 
tors localized upon each neighbourhood J\f x ; 

are two isometric operators localized upon 
each node x, and whose actions depend only on 
&mx(H x ). 



Moreover: 

• If the (TL X ) are all of finite dimensions, then the 
(K X ),D> and E are finite dimensional operators; 

• V U(® \q)) = (® |?», then \cf>) = ((g) \q)); 

• If the (TL X ) are all of infinite dimensions and 
U(®\q)) = (&)\q)), then we can choose to just 
have (®£>)((g) K X )(®E) = U\ip) where D and E 
are also unitary. 

Proof. [Encoding]. The action of E upon node x is just 
to add an ancilla, i.e. E\ip x ) = \q) <E) \ip x ). Hence if 
dim(W x ) is finite then E : H x — > H x ® H x and E is an 
isometry, whereas if TL X is of infinite countable dimension 
then we can use any bijection from N x N to N so that 
E : H x — > H x and E is unitary. 

[Product states]. Let us consider £Tt having the form 
of a product state, i.e. so that \ijj) — §Z)\ip x )- We will 
show that ((g) Z^X® K x )(® E)\ip) = U\ip), and then 
by linearity the result will be proved for entangled states 
also. This is because in general any state vector \<j>) can 
be written as a sum of | </>*), where each \(j) 1 ) is a product 
state \4> l ) = (g) \4> x ). Below we again use this form for 
10) =tf t (®l<7»- 

[Two tapes]. So E takes \ip) into ((g) |g))<8>(® \ipx))- Now 
since ((g) \q)) = UU^(§Z) \q)) = U\4>) we rewrite E\ip) as: 

i 

So initially our QLG has got two "tapes" , one which we 
call the "computed tape" holding state ?7((g \<t>)), and 
one which we call the "uncomputed tape" holding state 

(®hk». 

[Changing factors]. Now the idea is that the K x will let 
us pass pieces of the uncomputed tape to the computed 
tape. Namely we want K x E\ip) is equal to: 

<g) |0t))®(|^)<8 \4v)). 

V\{x} V\{x} 

Let us simply take K x = U Swap x W , meaning that we 
simply uncompute the computed tape, swap for \q), 
and then compute it back. Clearly this does the job but 
does seem wrong, because it looks as though we are act- 
ing over the entire graph and not just Af x . Yet this naive 
choice is actually the right one. Indeed since U is unitary 
causal with respect to T, then so is W with respect to T T , 
by Proposition^ And now since W is unitary causal over 
r T , U Swap x W must be localized upon M x , by virtue 
of Proposition [T] Note that the (K x ) commute with one 
another just because the (Swap x ) commute with one an- 
other and A i— > UAW is a morphism. 
[Decoding] . Of course we can reiterate this process until 
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we get 

£c/((g)l^>)®((g)l^>) 

i 

which is just U\ip) <g> \<j>). Now we just need to swap the 
computed and uncomputed tapes to get (|</>) ® U\tp). (In 
situations where t/^((g) \q)) = \<f>) is known and turns out 
to be a product state &)\<f>x}, then D can also locally 
undo the \</>) so as to get ((g) |g)) ® U\ip) = EU\tp). This 
is the case for instance in the standard situation when 
\q)) = ((g) \q))- If on top of that E was a unitary, 
D can also apply E^ and give back U\ip).) □ 

Corollary 1 (Circuit representation) 

Let U be a unitary causal linear operator with respect to 
a quantum labeled graph T = (V,£,7i). Then U can be 
expressed as a circuit of quantum operations each local- 
ized upon a neighbourhood M x , and having depth less or 
equal to deg(r) 2 + 2. 

Proof. By inspection of the proofs of Theorem [T] and 
using the following remarks. Since each K x is localized 
upon J\f x , many of them can be done in parallel, namely 
whenever the corresponding neighbourhoods do not in- 
tersect. The question of how much can be done in paral- 
lel, i.e. how many layers of circuit are necessary, is equiv- 
alent to the L(l, l)-labeling problem for graphs, namely 
we want to colour the graph so that no neighbours nor 
next-neighbours have the same colours. This is known to 
require at most deg(r) 2 colours Q. The plus two is for 
E and D. □ 
The study of unitary causal operators has older origins 
than the rise quantum information processing, for similar 
questions are clearly coming up in axiomatic/algebraic 
quantum field theories 0] — as argued also in the papers 
which treat the two-systems and three-systems cases of 
this theorem 0, [24|. The main difference in approach 
seems to be that AQFT caters looks at continuous time 
and space. The authors are not aware, however, of a 
result akin to Theorem [1] in AQFT, which would let us 
structure the dynamics of the system in such a meaning- 
ful, operational manner. But since our initial motivation 
was the study of QCA — let us see what this result has 
to say about them. 

V. CONSEQUENCES AND RELATED RESULTS 
IN QCA 

Cellular automata (CA) as introduced by Von Neu- 
mann psH ], consist of an array of identical cells, each of 
which may take one in a finite number of possible states. 
The whole array evolves in discrete time steps by iter- 
ating a function G. Moreover this global evolution G 
is shift-invariant (it acts everywhere in the same way) 
and local (information cannot be transmitted faster than 
some fixed number of cells per time step). Because this is 



a physics-like model of computation, Fcynman s ugg ested 
right from the birth of quantum computation [111 ] that 
one should look into quantizing this model (for two rea- 
sons: first because in CA computation occurs without 
extraneous control, hence this gets rid of a source of de- 
coherence; second because they are a good framework to 
study quantum simulation of a quantum system). Quan- 
tum cellular automata were always going to be unitary 
operators G over arrays of finite-dimensional systems, 
together with shift-invariance ("the laws of physics are 
everywhere the same") and causality ("there can be no 
instantaneous long-range communication" ) . At an infor- 
mal level this concept has been around for almost twenty 
years, unchallenged but yet somewhat impractical — in 
the sense that there was no proper axiomatization and 
nor a generic operational description of them. And so 
it was not really know what these things actually looked 
like. As Gruska puts it in one of the very first textbook 
on quantum computation [l^]: "A suitable definition of 
two- and more-dimensional quantum cellular automata 
is an untrivial matter." 

This situation has led to several competing definitions of 
QCA, each one attempting to tame the structure of the 
unitary operator G in its own manner. Let us briefly look 
at the three main approaches towards defining QCA — 
excluding by lack of space those works which are more 
concerned with quantum walks [Til . Il9| , quantum simu- 
lations [1, 01| , or implementation models [H, [27| ■ His- 
torically the first approach [l], [l^, Hil was recently shown 
to break causality [2j, and so it seems we must abandon 
this definition. The second approach [23[ is the axiomatic 
one, it provides a rigorous axiomatics for quantum cellu- 
lar automata, which we can rephrase in the vocabulary 
of this paper as follows: 

Definition 6 (QCA) 

A n-dimensional quantum cellular automaton (QCA) is 
an shift-invariant unitary causal operator G : TL — ► TL 
over a QLG Y = (V,£,TL) with: 

- V = Z™ i.e. the nodes form a grid; 

- £ = {x, x + z | x £ Z" A z £ {0, 1}"} i.e. radius half; 

- TL = (TLs) i-e. all cells are of a given finite dimension 
d. 

This paper [23j then derives a Block structure in a very 
general fashion. Unfortunately the proof is flawed in 
n-dimensions Q. Fortunately the result contained in 
this paper entails another Block representation of n- 
dimensional QCA: 

Theorem 2 (n-dimensional QCA) 

Let G be an n-dimensional QCA with alphabet S. Let E 
be an isometry from TLy: — > TLy, ® 7is such that E\ij) x ) = 
\q) ® l^x}- This mapping can be obviously extended to 
whole configurations, yielding a mapping E : Hcf ~^ 
Ji nTj i . Then there exists a n-dimensional QCA H on 
alphabet S 2 ; such that HE = EG. and H admits an 2 n - 
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layer block representation. Moreover H is of the form 

H = (®S)(Y[K X ) (2) 

where: 

• (K x ) is a collection of commuting unitary opera- 
tors all identical up to shift, each localized upon 
each neighbourhood Af x ; 

• S is the swap gate over We ® He , hence localized 
upon each node x. 

Proof. By inspection of the proof of Theorem [1] and us- 
ing the following remarks. At each x = . . . , i n ) step, 
K x = Ki lM ..i n is local to cells {ii,ii + l}x. . .x{i„, i n + l}, 
uncomputed and computed tapes alike. Namely, when- 
ever (ix, ■ . ■ , i n ) and . . . ,j n ) are such that for every 
k e {l,...,n}, \ik—jk\ > 1, then K^...^ and K jl ... Jn 
can be performed in parallel. So we can first apply 
simultaneously all the i n 's where the ifc's are even. 

Then, as each element x = (x\, . . . , x n ) can be written 
in a unique way as the sum of y with even coordinates 
and z S {0,1}", we need |{0, 1}™| = 2™ layers to 
apply all of the -Kii,...,t n 's. Moreover by shift-invariancc 
these ifj 1 ...j n 's are just shifted version of the same 
K, so that each layer is just tiling of the space by a 
finite-dimensional unitary K . □ 
(Notice that even the shift QCA can be implemented as a 
quantum circuit according to this theorem, but of course 
this is at the price of introducing ancilla. This question 
of the special role of the shift and its implementability 
is emphasized in Q and carefully analysed in [HI-) 
The third approach to define QCA is to give them an a 
priori hands-on, operational description of a particular 
form. Several works have followed this route for instance 
@ HHj H2, HH , but amongst them Q stands out at this 
stage as it just directly posits, after some interesting 
arguments, that their evolutions takes a form akin to 
the one described in given by Theorem [2j 
Here we have demonstrated that starting just from the 
axiomatic definition of QCA as in [23| and 0, one can 
derive a circuit-like structure for a QCA, thereby extend- 
ing the result of [23| to the n-dimensional case. We have 
also demonstrated that the constructive/operational 
definition of Q can be given a rigorous axiomatics. 
And by doing these two things we have shown that the 
definitions of [23[ and are actually equivalent up to 
ancillary cells. This clearly reinforces the feeling that 
the community has now got a well-axiomatized and yet 
concrete definition of n-dimensional QCA. 
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and struggle to obtain a practical representation for it. 
Often, however, these infinite dimensional Hilbcrt spaces 
arise from a position degree of freedom. By virtue of 
the principle according to which information travels at 
bounded speed, we can then think about "cutting space 
into different pieces" such that at each time step, the 
state of a piece depends solely on that of its neighbours. 
Whenever this happens Theorem [T] applies and lets you 
write U as (0 D f )(® if x )(® E), where D and E are 
local to each piece and K x is local to piece x and its 
neighbours. 

And so this is saying something very general which can 
be summarized by "Unitarity plus causality implies local- 
izability" . If a global evolution is locally implcmcntablc, 
this means we can focus on understanding local inter- 
actions between physical elements, and then the global 
evolution will just turn out to be a composition of them. 
And so this statement bridges a certain gap between gen- 
eral physical principles and the study of elementary inter- 
actions. Unfortunately there were two strongly limiting 
assumptions underlying this result, which we would be 
interested to lift. 

The first one is unitarity. A not so uncommon belief 
amongst theoretical physicists is that the universe being 
a closed system it should evolve unitarily. Nevertheless 
this is clearly an unpractical view — any everyday physi- 
cal system is an open system, noisy due to its interactions 
with the outside world, amongst which any measurement 
we may wish perform upon the system. Hence one of 
our most wanted open problem at the moment is to ex- 
tend Theorcm[T]to quantum operations. Clearly the way 
causality was axiomatized here will not entail localizabil- 
ity in an open systems setting (cf. PR-Boxes etc. 0]), 
so part of the challenge is to come up with a reinforced, 
yet intuitive notion of causality. Several cases of QCA 
under noise have been studied in [f| 0j| [l7| , which might 
perhaps be a guidance. 

The second one is discreteness. Clearly a point which is 
also rather open to discussion is whether we arc indeed 
allowed to divide up the universe into different "places" 
whose state depend only on that of the closest neighbour, 
in the sort of abrupt and discrete manner which we use 
here. Hence we would like to study the relationship be- 
tween QCA and continuous-space continuous-time mod- 
els, maybe building upon what has been done for quan- 
tum simulation from Quantum Lattice Gas Automata 

More generally it is our intention to understand the ex- 
tent in which "causality implies localizability" could be 
made into a general principle. Such a principle would 
impact theoretical physics, by providing operational de- 
scriptions of global evolutions in physics. But it would 
also impact theoretical computer science, as these op- 
erational descriptions become closer and closer to being 
computable descriptions of global evolutions in physics. 



There are many situations in physics where we want to 
study a unitary operator U over a large Hilbert space Ti, 
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